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We study decoherence of a three-qubit array coupled to substrate phonons. Assuming an initial 
three-qubit entangled state that would be decoherence-free for identical qubit positions, allows us 
to focus on non-Markovian effects of the inevitable spatial qubit separation. It turns out that the 
coherence is most affected when the qubits are regularly spaced. Moreover, we find that up to a 
constant scaling factor, two-qubit entanglement is not influenced by the presence of the third qubit, 
even though all qubits interact via the phonon field. 



I. INTRODUCTION 

A major obstable on the way towards a working quan- 
tum computer is decoherence: the interaction of the 
qubits with their environment reduces the indispens- 
able quantum coherence of the quantum states. Several 
strategies are pursued to beat decoherence. An active 
strategy is quantum error correction, which requires a 
redundant encoding of a logical qubit by several physi- 
cal qubits [3, H, [1] . Standard error correction protocols 
presuppose that all physical qubits couple individually 
to uncorrelated baths. A passive strategy is the use of 
decoherence-free subspaces (DFS) [3, [sTIE 0|- There, 
one logical qubit is encoded by several physical qubits in 
such a way that the logical qubit states do not couple to 
the environment. Ideal DFSs occur when physical qubits 
couple via a collective coordinate to a common bath. 

For solid-state qubits, the coupling to substrate 
phonons often is a relevant source of decoherence, in par- 
ticular for charge qubits in quantum dots Q. Whether 
these qubits experience correlated or uncorrelated noise 
depends on their distance in relation to the coherence 
length of the phonons, the sound velocity, the cutoff fre- 
quency and also on the dimensionality of the substrate. 
In Ref. @, this dependence has been worked out by 
studying pure dephasing of a two-qubit state with an ini- 
tial entanglement that is decoherence-free if both qubits 
couple to the environment at the same position [lOj. If 
the qubits are spatially separated, however, this behavior 
changes: the entanglement decays until the transit time 
of a sound wave from one qubit to the other is reached. 
However, if the qubits are embedded in a quasi-one- 
dimensional environment, the noise at the two positions 
eventually becomes sufficiently correlated to bring deco- 
herence to a standstill. In this way, a decoherence-poor 
subspace can emerge. Similar results can be found for 
the decoherence of entangled states of a regularly spaced 
qubit array [Til ] . 

This work is motivated by two questions: First, do im- 
perfections in the regular alignment of the qubits involve 
additional decoherence, or, put differently, how do irreg- 
ularities in the spatial extension of a qubit register influ- 
ence the collective decoherence properties? And second, 
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FIG. 1: Sketch of three qubits (green boxes) in a linear ar- 
rangement at positions x v , v = 1,2, 3, with the distances a + 8 
and a — 8. They interact via a coupling to the bosonic field 
(red line). 



is the entanglement of a qubit pair affected by its indirect 
interaction with a third qubit via the environment? 



II. QUBITS COUPLED TO A BOSONIC FIELD 

As sketched in Fig. [TJ we consider a linear arrange- 
ment of three qubits at positions x\ = 0, xi = a + 5, and 
X3 = 2a, i.e. the nearest-neighbor separations X12 = a + S 
and X23 = a — 6. To elaborate on the impact of spatially 
correlated noise we assume the qubit array to be em- 
bedded in a channel-like structure, as may be realized in 
carbon nanotubes or in linear ion traps. Thus, we treat 
the bosonic environment as effectively one-dimensional. 
The total Hamiltonian modelling this situation reads 
H = H q + H qh + H h , where H q = J2 v =i h£l v o vz l2 de- 
scribes three qubits with energy splittings h£t v with a vz 
being a Pauli matrix for qubit v = 1, 2, 3. Note that there 
is no direct interaction between the qubits. The bosonic 
field described by H\> = J^k ^k&l&fc consists of modes 
k with energies fuvk and the respective annihilation and 
creation operators bk and b\ . We assume a linear disper- 
sion relation ujk = c\k\ with sound velocity c. The transit 
time of a field distortion between the qubits v and v' is 
then t vv i = x vv i /c with aw — \x v — x u >\. Qubit v cou- 
ples linearly via the operator X v to the field, so that the 
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FIG. 2: Time evolution of the fidelity (|9} for the temperatures 6 = k^T /hui c = 10 4 (a) and 9 
coupling strength a = 0.001, and various displacements 8 of the middle qubit (see Fig. [TJ. 



10~ 3 (b), the qubit-field 



coupling Hamiltonian reads 

2 

tf qb =»X; X„£„, (1) 

with £„ = £ = ^ fc gke lkx " (bk + b_ k ) the bosonic field 
operator at the respective qubit position x v . We assume 
the microscopic coupling constants gk to be real-valued, 
isotropic, and identical for all qubits, i.e. gk v = gk and 
9-k = 9k- They determine the spectral properties of 
the bath and show up in the spectral density J(to) = 
J^kdk^i 10 ~ c ^)- Here we consider an Ohmic spectral 
density [13] 

J(Kj) = aue" B / Wt , (2) 

where the dimensionless parameter a denotes the overall 
coupling strength and lu c the cutoff frequency of the bath 
spectrum. 

The dynamics for the total density operator R of the 
qubits plus the environment is governed by the Liouville- 
von Neumann equation 

ihj t R(t)=[H qb (t),R(t)]. (3) 

The tilde denotes the interaction-picture representation 
with respect to H = H q + H h , i.e. A(t) = U^(t)AU (t), 
with time-evolution operator Uo(t) = exp{—iHot/h}. 
We assume that at time t = 0, the qubits can be pre- 
pared in a well-defined initial state, uncorrelated with 
the thermal bath. This constitutes an initial condition 
of the Feynman- Vernon type, where the total initial den- 
sity matrix i?(0) is a direct product of a qubit and bath 
density operator, i?(0) = p(0) <g> p^ 1 . The canonical en- 
semble of the bosons at temperature T is denoted by 
= exp(— Hb/kBT)/Z, with Z the partition function. 
We are interested in the reduced density matrix of the 
qubits p{t) — trb R(t) , where trb denotes the trace over 
the bath variables. 



III. DEPHASING AND ENTANGLEMENT 
DECAY 

In order to exemplify the impact of a spatial qubit 
separation on decoherence, we consider as the initial state 
the three-qubit entangled W state 

|W) = -^(|100> + |010> + |001>), (4) 

i.e. p(0) = |W / )(VF|, with the computational basis 
{|nin 2 n3)} where a vz \nin 2 n 3 ) = (-l) n "\nin 2 n 3 ) and 
n v = 0,1. Our motivation to focus on the initial state 
(U]) is twofold: First, W states play an important role 
in several protocols for quantum information processing 
HE [3 Ell, so that the entanglement dynamics after their 
preparation is relevant in itself. Second, the W state 
is special since it stays robust under collective dephas- 
ing , i.e. for vanishing qubit separations (x\ = x 2 = X3) 

Pure phase noise will also be assumed in the following, 
in which case the coupling operators in Eq. ([T]) become 
X v = <y vz . As a consequence, the interaction-picture 
qubit operators remain time- independent, X v (t) = X u . 
The exact time evolution of the reduced density operator 
can then be obtained, e.g. by a direct solution of the 
Liouville-von Neumann equation ^ [ll| . Amazingly, the 
exact result can even be obtained with an approximative 
time-local master equation approach already in second 
order of the qubit-field coupling a [l6| . It turns out that 
the density matrix elements in the basis {\n1n2n3)} at 
time t are proportional to their initial values. Thus, all 
density matrix elements that are initially zero remain 
zero, so that for the state \W), the dissipative quantum 
dynamics is restricted to the states 

|1) = |100), |2) = |010), |3> = ]001> , (5) 

i.e. at most nine out of 64 density matrix elements are 
nonvanishing. Initially they are all equal, i.e. p^y (0) = 
(j\p(0)\f) = 1/3, with = 1,2,3. They evolve as 

hi' (*) = \ exp{-A, y (t) + i[<f>j (t) - fa, (i)]} , (6) 
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FIG. 3: Final fidelity Fs(oo) as a function of the displacement 
S of the middle qubit, scaled to the value for equidistant qubit 
arrangement 5 = for various temperatures 9 = k&T /huj c and 
coupling strength a — 0.001. 



where the real part Ajj> (t) of the exponent accounts for 
the time-dependent amplitude damping of the matrix ele- 
ment; indirect interactions among qubits via the environ- 
ment give rise to a time-dependent frequency shift and 
the concomitant phase 4>j (t) — (t) with cf>j (t) = <pj (t) — 
Ailjt. Here Afij is a static frequency renormalization 
flj —> Qj + Ailj for qubit j and ipj(t) describes its onset; 
cf. Ref. UJ. Henceforth we work in the interaction pic- 
ture with respect to the renormalized energies. With the 
scaled temperature 9 — k^,T /huj c and the scaled times 
t = Lo c t and t vv i — uj c t vv i, the density matrix elements 
become £,y(r) = exp{-A jy (r) + i[ipj(r) - ^'(t)]}/3 
[111 ] with the phases 



<Pi(r) = ~~ E (-1) S '" +S -' E arctan t T ± ' ( 7 ) 



and the amplitude damping exp{— Ajj/(r)} 
f{T,Tjj')/f(T,0), where 



f(r,r') = 



|r(^[l-ir'])| 



16a 



\T 2 (6[1 + i{ T - r')])r 2 (6»[l + *(t - t> 



1 4a 



(8) 



|1 + r 2 [l - i T ']-2|4« 



As expected for pure dephasing, populations are pre- 
served: Pjj{t) — pjj(0). This implies that neither the 
qubits nor the total system will reach thermal equilib- 
rium. Nevertheless decoherence does occur since relative 
phases between eigenstates will be randomized so that 
off-diagonal density matrix elements decay. 

Only in the absence of the environment (a = 0), the 
qubits remain in the W state (dJ). As a measure for the 
deviation from this "ideal" output state p(0) = |T4 7 )(H 7 |, 
we employ the fidelity F(t) = tr{p(t)p(0)} [lj], which in 
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(9) 
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The index S refers to the displacement of the middle qubit 
v = 2. The time evolution of the fidelity is shown in 
Figs. [2Kb for two different temperatures. Clearly, the 
fidelity decay is slowed down whenever a transit time t = 
tjji is reached. At time t — tis, when the field has also 
enabled communication between the two outer qubits 1 
and 3, decoherence even comes to a standstill! Note that 
other initial states may lead to complete dephasing [l7l ]. 
The fidelity saturates to a finite value Fg(t > t 13 ) = 
Fs (oo) which is larger the lower the temperature. For 
a fixed temperature, the stable fidelity increases if the 
middle qubit is displaced away from the central position 
and is maximal for 6 = ±a, i.e. if qubit 2 becomes co- 
located with qubit 1 or 3. The fidelity gain as a function 
of the asymmetry 5 /a is shown in Fig. [3] 

The explanation of this behavior of Fs (oo) follows from 
Eq. @. The individual coherences p VV ' decay approxi- 
mately exponentially and with the same decay rate [ill ], 
but stop decaying at different times \x u — x v i\/c. Hence 
for the fidelity it pays off to displace the middle qubit, 
thereby stopping the decay of P23 at an earlier time 
(a — S)/c at the expense of stopping the decay of P12 
only at time (a + S) / c. 

For the extreme cases 8 — and S — a, the fi- 
delity decreases monotonously with temperature, but for 
< 5 < a an optimal temperature exists for which the 
fidelity gain is maximal. This behavior is related to the 
fact that for all temperatures pi2(i) = P23(t) when 5 = 0, 
and pi3(t) = P2a(t) when S = a. Only in the general 
asymmetric configuration do all three coherences stop de- 
caying at different times. 

The initial W state is special in the sense that it ex- 
hibits bipartite entanglement between any qubit pair. 
But does the middle qubit affect the entanglement de- 
cay of the outer ones in any way? The reduced density 
matrix of qubits 1 and 3 is 



fr2/3(t) = o (]00><00j + |01)(01| + |10)(10|) 



■Pi3(*)|10><01| +p 3 i(t)|01><10] 



(10) 



It depends on time only through the coherence pi3(t) = 
pjji(i) and thus, according to Eq. ([6]), only on the transit 
time ti3 = 2a/c between the outer qubits. The entangle- 
ment of the outer qubits is therefore independent of the 
displacement S of the middle qubit! Indeed, for their con- 
currence we find C[tr 2 p(t)] = 2|pi 3 (i)|. The same 
dynamics is found for the two-qubit concurrence when 
only the outer two qubits had been present, with the 
two-qubit W state (i.e. the robust Bell state) as their 
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initial state @ . Surprisingly, the only effect of the pres- 
ence of the middle qubit is a rescaling of the concurrence 
of the outer two by a factor 2/3. 

IV. CONCLUSIONS 

We have studied the pure dephasing of three spatially 
separated qubits in an Ohmic environment. Dephasing 
can be incomplete when starting in a W state, not only 
for symmetrically spaced qubits. For fixed separation of 



the outer two qubits, the final fidelity is even larger the 
more "clustered" the qubits are. Surprisingly, the middle 
qubit does not affect the dynamics of the concurrence of 
the outer two. 
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